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ABSTRACT 

The observational signatures of many of the most fundamental physical attributes of galaxies, 
including the mass of the interstellar gas in both atomic and molecular form, the thermody- 
namical properties of this gas, as well as the ambient magnetic fields and cosmic rays depend 
on the strength and colours of the ambient radiation fields, which however are generally not 
directly observable in star-forming galaxies. Here we provide and describe a library of diffuse 
stellar radiation fields in spiral galaxies derived using calculations of the transfer of stellar 
radiation from the main morphological components - disks, thin disks, and bulges - through 
the dusty interstellar medium. These radiation fields are self-consistent with the solutions for 
the integrated panchromatic spectral energy distributions (SEDs) previously presented using 
the same model. Because of this, observables calculated from the radiation fields, such as 
gamma-ray or radio emission, can be self-consistently combined with the solutions for the 
ultraviolet/optical/submillimeter SEDs, thus expanding the range of applicability of the ra- 
diation transfer model to a broader range of wavelengths and physical quantities. We find 
strongly contrasting solutions for the spatial distribution of the radiation fields for disks, thin 
disks and bulges, which are also strongly dependent on dust opacity r B , and wavelength. We 
give an analytic formula for the optically thin case of an infinitely thin disk, which we show 
well approximates the solutions obtained for the thin disk at longer optical wavelengths. For 
bulges we find a strong dependence of the radiation fields on Sersic index. Consequently, we 
show that dust heating in the central regions of galaxies with different bulge Sersic indices 
can vary by more than a factor of 10. For resolved galaxies with high Sersic index bulges, 
this results in the appearance of a prominent mid-infrared source in the central regions of 
disks, powered by optical photons. All the solutions from the library of radiation fields which 
contains 630 two-dimentional spatial grids (in cylindrical coordinates) of energy densities are 
made available in electronic format. 

Key words: radiative transfer - dust, extinction - galaxies: ISM - galaxies: spiral - galaxies: 
stellar content 



1 INTRODUCTION 

Knowledge of radiation fields (RF) in galaxies is of prime impor- 
tance to many branches of astrophysics, embracing low- and high- 
energy processes in the interstellar and intergalactic medium, with 
applications to cosmology, galactic astrophysics, astro-chemistry 
and astro-biology. For example ultraviolet (UV) radiation fields 
are one of the most important heating agents for interstellar gas, 
proceeding either through photoionisation in the presence of ion- 
ising UV radiation, or through the photoelectric effect in the neu- 
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tral ISM. Thus, the main cooling lines of the interstellar medium, 
in particular the [CII] line, are produced as a result of the photo- 
electric effect on dust grains, primary driven by non-ionising UV 
photons from the ambient radiation fields interacting with poly- 
cyclic aromatic hydrocarbon (PAH) molecules and dust grains, both 
in star-formation regions (e.g. Crawford et al. 1985, Helens & 
Hollenbach 1985, Wolfire et al. 1989, Hollenbach 1991, Israel & 
Maloney 2011) and in the diffuse ISM (e.g. Madden et al. 1993, 
Bennett et al. 1994, Pierini et al. 1999, Pierini et al. 2001, Rubin 
et al. 2009, Lebouteiller et al. 2012). Diffuse RFs impinging on 
molecular clouds in the interstellar medium also drive the exita- 
tion of molecules in the outer layers of the clouds, determining the 
rotational-line emission from rare tracers such as CO molecules, 
from which the bulk of the mass in the form of molecular hydrogen 
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can be determined. A quantitative knowledge of the amplitude and 
spectrum of the RFs is also a necessary input for hydro-dynamical 
simulations of the turbulent structure of the interstellar medium 
(e.g. de Avillez & Breitschwert 2007). 

In high-energy astrophysics it is particularly important to 
know the number density of all photons, including those of low 
energies where the number density is high, because such photons 
can be transformed into gamma-rays through the inverse Comp- 
ton scattering (ICS) with relativistic electrons. Indeed, inverse 
Compton scattering provides one of the principal y-ray produc- 
tion mechanisms (Jones 1968, Blumenthal & Gould 1970, Aha- 
ronian & Atoyan 1981, Narginen & Putanen 1993, Brunetti 2000, 
Sazonov & Sunyaev 2000). In a variety of astrophysics environ- 
ments, from very compact objects like pulsars and active galactic 
nuclei to extended sources like supernova remnants and clusters of 
galaxies, low-energy photons are efficiently boosted to high ener- 
gies through this mechanism (Aharonian & Ambartsumyan 1985, 
Zdziarski et al. 1989, Dermer & Atoyan 2002, Moderski et al. 
2005, Khangulyan & Aharonian 2005, Kusunose & Takahara 2005, 
Stawarz et al. 2006, 2010). In galaxies the UV/optical and far- 
infrared (FIR) interstellar radiation fields can be an important chan- 
nel for the production of inverse Compton y-rays. Thus, knowl- 
edge of the interstellar RFs is needed to measure the distribution 
of relativistic electrons in these objects, as well as to isolate other 
gamma-ray emissions, such as those produced in the decay of pi- 
ons generated in the interaction of relativistic protons with interstel- 
lar gas. Similarly, because relativistic electrons also produce radio 
synchrotron emission, knowledge of the RF in galaxies is neces- 
sary to elucidate the strength of interstellar magnetic fields from a 
combination of radio and gamma-ray observations. Radiation fields 
also affect the propagation of gamma-rays through the intergalactic 
medium (e.g. Stecker et al. 1992) as well as even in the interstel- 
lar medium (Moskalenko et al. 2006) through pair production of 
gamma-ray photons on the background radiation fields. 

In summary it is clear that the observational signatures of 
many of the most fundamental physical attributes of galaxies, in- 
cluding the mass of interstellar gas, in both atomic and molecular 
form, the thermodynamical properties of this gas, as well as the am- 
bient magnetic fields and cosmic rays in the form of both protons 
and nuclei and electrons, depend on the radiation fields in galax- 
ies. It would be relatively straightforward to determine radiation 
fields by direct observations of the sources of photons, which, in 
most cases, are predominantly stars, if galaxies were optically thin. 
However, galaxies - and in particular star-forming galaxies - con- 
tain dust, which, because it absorbs and scatters photons, partially 
or wholly prevents a direct measurement of the spatial and spectral 
distribution of the sources of the RFs. The propagation of light de- 
pends in a complex way on the relative distribution of stellar emis- 
sivity and dust opacity, with structures ranging in scale from par- 
sees to kiloparsecs. Fortunately, although dust is the cause of this 
problem, it is also part of the solution, because the absorbed light is 
re-radiated in the mid-infrared (MIR)/FIR/submillimeter (submm). 
In particular grains which are sufficiently large to be in equilibrium 
with the radiation fields directly probe the strength of these radi- 
ation fields, acting as a strong constraint on the radiation transfer 
problem. In addition, the radiation emitted by grains small enough 
to be impulsively heated by individual photons strongly depends 
on the colour of these photons, placing particular constraints on the 
UV component of the interstellar radiation fields, providing one has 
independent knowledge of the size distributions of the constituent 
populations of grains. 

Thus, the key to determining radiation fields in galaxies is 



to invert the observed broad-band spectral energy distribution in- 
corporating both measurements of direct light in the ultraviolet 
(UV)/optical, as well as the dust re-radiated light in the Far-infrared 
(FIR)/submillimeter (submm), self-consistently taking into account 
constraints on the spatial distribution of both stars and dust as well 
as on the optical properties of the grains. This is the only way in 
which the exact contribution of optical and UV photons in heat- 
ing dust grains can be traced, allowing for a proper determination 
of the contribution of old and young stellar populations in heat- 
ing the dust (Popescu & Tuffs 2010), and therefore allowing for a 
consistent estimation of the star-formation rates in galaxies, and, 
ultimately for a robust determination of their star-formation his- 
tory. Self-consistent modelling of the panchromatic spectral energy 
distributions (SEDs) of galaxies (e.g. Silva et al. 1998, Popescu et 
al. 2000, Misiriotis et al. 2001, Popescu et al. 2004, Siebenmorgen 
& Kriigel 2007, Bianchi 2008, Baes et al. 2010, Baes et al. 2011, 
Popescu et al. 201 1, MacLachlan et al. 201 1, Schechtman-Rook et 
al. 2012, de Looze et al. 2012a,b) have been primarily motivated by 
the need to extract intrinsic parameters of the stars and dust, most 
notably the aforementioned SFRs, SF histories and dust masses. 

In previous papers (Popescu et al. 2000, Tuffs et al. 2004, 
Popescu et al. 2011) we provided self-consistent solutions for 
modelling the integrated panchromatic SEDs of star-forming disk 
galaxies. Here we give and describe the solutions for the radiation 
fields, implicit to the SEDs presented in our previous work. The so- 
lutions are given in the form of a library of stellar radiation fields for 
the main morphological components of our model of spiral galax- 
ies: disks, thin disks and bulges. These radiation fields are identical 
to those used when calculating the library of panchromatic SEDs 
presented in Popescu et al. (201 1). Because of this observables cal- 
culated from the radiation fields, such as gamma-ray or radio emis- 
sion, can be self-consistently combined with the solutions for the 
UV/optical/submm SEDs found using Popescu et al. (2011), thus 
expanding the range of applicability of the radiation transfer model 
to a broader range of wavelengths and physical quantities. Accord- 
ingly, the library of radiation fields is calculated for different values 
of central face-on dust opacity, as defined in Popescu et al. (201 1). 
The separate tabulation of the radiation fields for disks, thin disks 
and bulges is possible because all the stellar components are at- 
tenuated by the same distribution of dust, and therefore the diffuse 
radiation fields seen by each grain can be considered to be a sum 
of the radiation fields produced by each stellar component. Thus, 
the radiation fields can be combined according to other model pa- 
rameters from Popescu et al. (2011): the bulge-to-disk ratio (B/D), 
the 5 FR determining the luminosity of the young stellar population 
(thin disk), and the parameter expressing the normalised luminosity 
of the old stellar population in the disk, old, to derive combined ra- 
diation fields of spiral galaxies. This concept, which relies on stellar 
light being an additive quantity, was originally introduced in Tuffs 
et al. (2004) to describe the attenuation of stellar light. Here we 
only show radiation fields in direct stellar light (UV/optical/NIR), 
as these are the ones that can be separately calculated and have ad- 
ditive properties. Radiation fields in dust re-emission will be the 
object of a future paper, as these cannot be separated according to 
the different stellar components heating the dust, nor have they ad- 
ditive properties. This is because there is no linear transformation 
between dust emission and absorption. For all these reasons, the 
radiation fields in dust re-emission can only be presented as com- 
bined radiation fields, characteristic of a galaxy as a whole. 

This paper is organised as follows. The general characteristics 
of the calculations are presented in Sect. [2] In Sect. [3] we describe 
the spatial variation of the radiation fields, for each morphologi- 
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cal component. In Sect.[4]we analyse the variation of the radiation 
fields with face-on dust opacity. For the specific case of bulge radia- 
tion fields we also describe in Sect.[5]the change in energy densities 
when changing the Sersic index of the volume stellar emissivity. In 
Sect.[6]we present the library of radiation fields and show how the 
energy densities can be scaled according to the size and luminosity 
of a given galaxy. In Sect. [8] we give our summary. 



2 SIMULATIONS 

The simulated radiation fields are those used to generate the large 
library of dust- and PAH-emission SEDs presented in Popescu et 
al. (2011). This library is self-consistently calculated with the cor- 
responding library of dust attenuations first presented in Tuffs et 
al. (2004) and then in updated form in Popescu et al. (2011). The 
details of the calculations leading to the production of the the radi- 
ation fields are described at length in Popescu et al. (2011). Here 
we only briefly mention their main characteristics. All calculations 
were made using a modified version of the ray-tracing radiative 
transfer code of (Kylafis & Bahcall (1987)), which treats scatter- 
ing utilising the method of scattered intensities (Kylafis & Xilouris 
2005), and the dust model from ( TWeingartner & Draine ( 2001)) 
and ( |Draine & Li (2007)) , incorporating a mixture of silicates, 
graphites and PAH molecules. 

The radiation fields are calculated separately for the old stel- 
lar disk ("the disk") , the bulge and the young stellar disk ("the thin 
disk) of a galaxy, all seen through a common distribution of dust. 
The geometrical model of Popescu et al. (2011) consists of both 
of a large-scale distribution of diffuse dust and stars, as well as of 
a clumpy component physically associated with the star-forming 
complexes. For the puipose of this study only the large-scale dis- 
tribution of diffuse dust is considered, as it is this that affects the 
large-scale distribution of UV/optical light determining the values 
of the diffuse radiation fields. 

The intrinsic volume stellar distributions are described by ex- 
ponential functions in both radial and vertical direction for the two 
stellar disks and by deprojected de Vaucouleurs or Sersic functions 
for the bulges. The diffuse dust is distributed in two disks associated 
with the old and young stellar populations, which we refer to as the 
thick and thin dust disk, respectively, on account of the different 
thickness of these structures. These dust distributions are described 
by double (radial and vertical) exponential functions. The thick dust 
disk is less opaque than the thin dust disk, with the relative face-on 
central dust opacity being fixed to a value of 0.387. 

The length parameters of the model describing the volume 
emissivity for stars and dust - scale-lengths, scale-heights and ef- 
fective radii - are listed in Table 1 in Tuffs et al. (2004). The rele- 
vant information for this work is that the old stellar disk component 
has a scale-length that decreases with increasing optical/NIR wave- 
length, as given in Table 2 in Tuffs et al. (2004), while the scale- 
height remains constant over this wavelength range. Similarly, the 
effective radius of the bulge does not vary with optical/NIR wave- 
length. The bulge is an oblate ellipsoid with an axial ratio (thick- 
ness) of 0.6. The young stellar disk has a much smaller scale-height 
than the older stellar disk (by a factor of 4.6), while its scale-length 
is constant over wavelength and is equal to that of the old stellar 
disk in the B band. 

The scale-length of the thick dust disk associated with the old 
stellar population is larger (by a factor of 1 .4) than that of the cor- 
responding stellar disk, while its scale-height is smaller (by a fac- 
tor of 1.5) than the scale-height of the old stellar disk. In contrast, 



the young stellar disk spatially coincides with the associated thin 
dust disk (same scale-heights and -lengths). As we will see in the 
next sections, the different vertical structure of the two dust disks 
strongly influences the solution for the radiation fields. The stellar 
emissivities of the disks are truncated at a radius of five exponen- 
tial scale-lengths while bulges are truncated at ten effective radii. 
The geometrical parameters of the model have either been empir- 
ically constrained to follow the median trends found by Xilouris 
et al. (1999) for a sample of well resolved edge-on spiral galax- 
ies, including NGC891, or have been fixed from physical consid- 
erations. The physical interpretation of this model is described at 
length in Tuffs et al. (2004) and Popescu et al. (201 1). A schematic 
representation of the geometrical model can be found in Fig. 1 in 
l |Popescu etal. (2011)| >. 

The radiation fields are sampled on an irregular grid with sam- 
pling intervals in radial direction ranging from 50 pc in the centre 
to up to 2 kpc in the outer disk and with sampling intervals in ver- 
tical direction ranging from 50 pc in the plane to up to 500 pc in 
the outer halo. This scheme was necessary to properly sample the 
radiation fields in the more optically thick central regions, where 
the dust associated with the thin stellar disk in conjunction with the 
strong positional dependence of the emissivity of the young stellar 
population (in z) and the inner emissivity cusp of the bulge (in both 
R and z) can lead to quite small scale structures even in the diffuse 
radiation fields. We made some further tests with higher sampling, 
ranging from 25 pc to up to 1 kpc in radial direction and from 25 pc 
to up to 250 pc in vertical direction. This made a less than 1% differ- 
ence in the results with a large increase in the computational time. 
We have therefore concluded that the optimum sampling of the ra- 
diation fields is in the AR = 50 pc - 2 kpc and Az = 50 pc - 500 pc 
scheme. 

All the calculated radiation fields are made available in elec- 
tronic format. Their main characteristics concerning the spatial 
variation, as well as the variation with dust opacity, are described 
in the next two sections. 



3 SPATIAL VARIATION OF RADIATION FIELDS 

Here we describe the spatial variation of the radiation fields of each 
individual stellar component of our model. Since, as mentioned be- 
fore, the solutions for each component are additive, solutions for 
composite structures can be found by weighted addition, according 
to the intrinsic stellar luminosity of each component. All examples 
are for the case = 3.5. This value of dust opacity corresponds 
to the solution we found for NGC 891 (Popescu et al. 2011) and 
for statistical samples of galaxies (Driver et al. 2007, Grootes et al. 
2013). 

3.1 The thin disk 

3.1.1 The radial profiles 

The radial variation of the RF of the thin disk changes drastically 
between an optically thick and optically thin case. To illustrate this, 
we show in Fig. Q] examples of RFs calculated at different wave- 
lengths. We first describe the behavior of the RF in the plane of the 
disk (z = 0). 

i) The optically thin case; z = 0. This arises at longer wave- 
lengths, in the J and K bands, where the profile shows a monotonic 
decrease with increasing radial distance. For comparison, the so- 
lution with no dust at all is overplotted with a red line in the panel 
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R [kpc] 

Figure 1. Radial profiles of radiation fields from the thin disk, for = 3.5. The solid line represents the profiles in the plane (z = Opc) while the other 
lines (see legend) represent profiles at different vertical distance from the plane (z =50, 200, 1000 and 2000 pc). The different panels show profiles at different 
wavebands, from the UV to the NIR. The profiles have been scaled such that the maximum intensity of the radiation fields of the model galaxy at different 
wavelenths reaches the value of 7.42 X 10 38 erg/pc 3 /A, which is the maximum intensity at 912 A. The profile overplotted in red in the bottom right panel (the 
K-band case) represents the dustless solution. 



corresponding to the K-band. One can see that the radial profile cal- 
culated for the longer optical wavelengths and for = 3.5 tends to- 
wards the completely optically thin solution, which shows an even 
steeper decrease with increasing radial distance. This is because for 
a thin disk geometry the component of energy density arising from 
sources within galactocentric radius R is well approximated by an 
inverse-square law for a point source situated at R = and with lu- 
minosity equal to the total enclosed luminosity within R. We show 
that the exact analytic formula of the infinitely thin disk is given by 
the sum of two terms: 



u(R) = m(R) + u 2 (R) 
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where a is the luminosity surface density. In our case a decreases 
exponentially with R, with an exponential scale-length h s . The first 
term in Eq.[T]falls off more slowly than an inverse square law for 
a point source situated at R = and with a luminosity equal to the 
whole luminosity of the galaxy. The second term in Eq.Q~]gives the 
contribution from sources beyond R. This term gives a more slowly 
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R [kpc] 

Figure 2. Comparison between a radial profile of radiation fields from the 
thin disk with no dust calculated with our radiative transfer code and the 
one given by the analytic formula from Eq. [T]for an infinitely thin disk. In 
both cases the calculation is for the K band. 



varying component for the RFs than the first term, so also neces- 
sarily falls off more slowly than an inverse square law for a point 
source of the same total luminosity. Beyond the radius enclosing 
half of the total luminosity, there will be a transition towards an in- 
verse square law. To compare the solution of the analytic formula 
for the infinitely thin disk with the real case of a thin disk (small 
but non-zero thickness) with no dust we show in Fig.[2]a compari- 
son between the radial profiles predicted by the two cases. One can 
indeed see that the analytic formula is a good approximation of the 
optical thin case. The small differences between the two profiles 
which appear at low radii are due to the thickness of the thin disk. 

ii) The optically thick case; z = 0. This arises at shorter UV 
wavelengths, where one can see a flat profile in the inner region, 
followed by a monotonic decrease with increasing radial distance. 
The inner flat profile is a direct consequence of the optically thick 
character of the solution. Thus, in an optically thick case the ra- 
diation seen at each point is proportional to the stellar emissivity 
within the horizon sphere (the surface for which r = 1 from the 
point considered) divided by the dust opacity: 



where rj is the stellar emissivity and k is the extinction coefficient 
(defined as the probability of extinction of a photon per unit length). 
The horizon sphere is then the sphere with radius equal to 1 /k. If 
both the stellar emissivity and the dust opacity have the same spa- 
tial variation within the horizon sphere (which is the case for the 
thin disk), then the RF will remain constant. As the radial distance 
increase, the horizon surface starts to become a prolate ellipsoid, 
since light can propagate to larger distances vertically than in the 
radial direction before reaching the r = 1 surface (the scale-height 
of the dust is much shorter than its scalelength, by a factor 0.016 in 
our model). Because the height of this ellipsoid becomes larger than 
the scale-height of the thin stellar disk and increases progressively 
with increasing radius, the ratio of stellar emissivity to dust opacity 
(in the vertical direction) will start to decrease with radius, because 
dust opacity from the thick dust disk becomes significant. Thus, ac- 
cording to Eq. [8] RF will start to fall off from the maximum value 
reached in the (completely) optically thick case. Eventually the RFs 
will tend towards an inverse square law decrease, which sets in at 



radii beyond which all outwards lines of sight (in the plane as well 
as perpendicular to the plane) become optically thin. 

iii) Intermediate cases; z = 0. At longer UV/shorter optical 
wavelengths, the profiles show a mixed behaviour, between the op- 
tically thin and optically thick limits. Essentially they resemble the 
profiles from the optically thick cases, with the difference that the 
flat part in the inner radius is missing. This is due to the fact that 
already in the inner region the horizon surface is given by a prolate 
ellipsoid and not by a sphere. Because the height of this ellipsoid 
is always significantly larger than the scale-height of the thin stel- 
lar disk and increases progressively with increasing radius, there is 
a corresponding decrease in the ratio of stellar emissivity to dust 
opacity (in vertical direction), resulting in the monotonic decrease 
of the RFs. 

Secondly, we describe the behaviour of the radial profiles at 
various vertical distances from the plane. 

i) The optically thin case; high z. As before, this arises at 
longer optical wavelengths. The profiles show a similar behavior 
to the radial profiles in the plane of the disk, but with a shallower 
decrease. This happens because the fractional change in the dis- 
tances to all emitting elements is smaller from an outside view than 
from an inside view (within the disk). 

ii) The optically thick case; high z. The profiles at large ver- 
tical distance are rather flat, as seen in Fig. Q] in the UV range at 
z = 1 - 2 kpc. This qualitative change in the solution for high z 
is because the thick dust disk, rather than the thin dust disk, con- 
trols the emission escaping from the thin disk. Even from a vantage 
point 2 kpc above the disk (the maximum height calculated here), 
the thin stellar disk is viewed through an optically thick layer of 
dust in the thick disk extending over almost 2 n sterad. The level 
of the energy densities of the RF is determined by the ratio of thin 
disk emissivity to thick disk k at the depth of the emitting surface 
from which the escaping photons are emitted. Because of the much 
steeper z dependence of the thin disk emissivity compared to the 
thick disk dust opacity, this depth is well below the t = 1 surface 
(as seen from above)- ie the effective optical depth of this surface is 
> 1 . This in turn means that the emitting surface has a larger extent 
than the extent of the r = 1 layer as seen from z = 0, explaining 
the much flatter form of the RF curve with radius for the z = 2 kpc 
case compared to the z = kpc case. 

3.1.2 The vertical profiles 

The trends in the vertical profiles can be easily explained following 
the trends already discussed for the radial profiles. Thus, at shorter 
wavelengths, where the disks are optically thick, the profiles in the 
centre (R = 0) show a decrease in the energy density with increas- 
ing z, until a certain distance above the disk, from where the RF 
remain fairly constant (see Fig. [3}. As discussed for the radial pro- 
files, in an optically thick case, the RF vary according to Eq. [8] as 
long as the vantage point is within the dust layer. However, unlike 
the radial variation, in the vertical direction the RFs are affected 
by the two dust disk structure, where the disk opacity has a dif- 
ferent variation from that of the stellar emissivity. The RFs will 
therefore decrease with increasing vertical distance from the disk, 
and this behavior continues also for vertical profiles at larger radii, 
since the disk remains optically thick for a large range of radial dis- 
tances. Moving still further out in z, it becomes possible for light 
to freely escape upwards from an emitting surface (below) which 
marks the transition (in the vertical direction) from optically thick 
to optically thin. At a vantage point more than 500 pc up, the contri- 
bution to the RFs from this surface dominates the in situ contribu- 
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Figure 3. Vertical profiles of radiation fields from the thin disk, for = 3.5. The solid line represents the profiles going through the centre (R = Opc) while 
the other lines (see legend) represent profiles at different radial distances from the plane (R =3, 8, 15 and 24 kpc). The different panels show profiles at different 
wavebands, from the UV to the NIR. The profiles at different wavebands have been scaled as in Fig. [T]The profile overplotted with the red line in the bottom 
right panel (the K-band case) represents the dustless solution. 



tion from the thin stellar disk, whose emissivity is tiny at more than 
a few scale heights. Because at short wavelengths the thin disk is 
optically thick (face-on) over a large range of galactocentric radii, 
the emitting surface extends over almost 2/r sterad. Thus, a vantage 
point at larger vertical distances above the disk will essentially see 
a constant surface-brightness until large radii, causing the RFs to 
remain constant higher-up above the disk. 

The trend of decreasing the RF with increasing z followed by a 
flat trend with z is gradually changing from shorter to longer wave- 
lengths, as the disk becomes more optically thin. Eventually in the 
optically thin case there is only a continuous decrease of the RF 
with increasing z. For comparison, the solution with no dust at all 
is overplotted with a red line in the panel corresponding to the K- 
band. One can see that the vertical profile calculated for the longer 



optical wavelengths and for = 3.5 tends towards the completely 
optically thin solution, though the latter shows a much steeper de- 
crease with increasing vertical distance. 

Coming back to the optical thick cases, but looking at the pro- 
files at larger radii (Fig. [3} one can see a similar behavior to that in 
the centre of the galaxy. However, at larger vertical distances above 
the disk, where the energy densities become constant, the overall 
level of the constant is higher than that of the RF in the centre of 
the galaxy. So the profiles do not show a monotonic behavior in 
this respect. This is because, although the vantage point above the 
disk will also see a constant surface brightness, the level of this is 
mainly determined by the thin (small scale-height) dust disk, since 
the thick dust disk (with the larger scale-height) becomes more op- 
tically thin. In other words at smaller radii and high z the radiation 
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Figure 4. Radial profiles of radiation fields from the disk, for = 3.5. The solid line represents the profiles in the plane (z = Opc) while the other lines 
(see legend) represent profiles at different vertical distance from the plane (z =50, 200, 1000 and 2000 pc). The different panels show profiles at different 
wavebands, from the B band to the NIR. The profiles have been scaled such that the maximum intensity of the radiation fields if the model galaxy at different 
wavelenths reaches the value of 18.54 X 10 38 erg/pc 3 /A, which is the maximum intensity in the B-band. The profile overplotted in red in the bottom right panel 
(the K-band case) represents the dustless solution. 



fields are proportional to 

rim 



u(R) ~ 



Kl(.R) + K 2 (R) 



(7) 



while at intermediate radii and high z the radiation fields are pro- 
portional to 



u(R) 



K 2 (R) 



(8) 



where /c, and k 2 are the opacities of the thick and thinn dust disks, 
respectively. As a consequence, the surface brightness increases at 
a intermediate radii, since the layer of constant surface brightness 
moves towards the plane of the disk. There will be a radius at which 
a maximum in the value of the energy density of the RF at larger z 
is attained, after which this will start to slowly decrease again. This 
causes the non-monotonic behaviour of the radiation fields above 
the disk with increasing radius. At longer wavelength this never 
happens, because the thin dust disk becomes transparent even in 
the inner radii. 



3.2 The Disk 

In our model we consider the disk to be a source of only opti- 
cal/NIR photons from an older stellar population, with stars old 
enough to have been dynamically heated to a thicker disk configu- 
ration. 

At shorter optical wavelengths the radial profiles at z = are 
rather flat (see Fig. [4}. This behaviour is quite different from that 
of the thin disk at the same wavelengths (see Fig. [TJ and appears 
because the stellar disk has a larger scale-height than that of the thin 
stellar disk. Thus, as in the thin disk case, at small radii and at z = 
the disk is optically thick in the radial direction and more optically 
thin in the vertical direction. Since the distance to the r = 1 surface 



in vertical directions is larger than in radial directions, the horizon 
surface is given by a prolate ellipsoid. Although the height of this 
ellipsoid continuously increases with radius (as in the case of the 
thin disk), this still remains smaller than the scale-height of the 
stellar disk. So, unlike the thin disk case, the ratio between stellar 
emissivity to dust opacity will remain fairly constant within the 
galactocentric radius for which the height of the ellipsoid is smaller 
than the scale-height of the stellar disk. Therefore, radial profiles 
for the radiation fields is flatter in the inner disk. 

At larger radii though, the disk become optically thin, consis- 
tent with a radial profile of RF decreasing with increasing radius. 
At longer optical/NIR wavelengths, the disk becomes optically thin 
even in the inner part, therefore the RF will show a continuous de- 
crease with increasing radius. For comparison, the solution with no 
dust at all is overplotted with a red line in the panel for the K-band. 
One can see that the radial profile calculated for the longer optical 
wavelengths and for = 3.5 is very well described by the com- 
pletely optically thin solution. 

The behavior of the RFs above the disk can be better seen in 
the vertical profiles (see Fig. [5}- The dominant feature of these pro- 
files is the existence of a maximum in the energy density of the 
RFs at an intermediate distance above the plane. Thus, the RFs in- 
crease from z = until a certain height, and then they decrease 
again. This is because the stellar emissivity in the disk has a larger 
scaleheight than the scaleheight of both dust disks. So in the opti- 
cally thick regime the energy density of the RFs goes up until the 
solution for the optically thin case is reached. The height for which 
the maximum is obtained decreases continuously with increasing 
optical wavelength (due to the decrease in the optical thickness of 
the disk). In the K-band the maximum is obtained at very small dis- 
tances above the plane, but it never reaches the z = plane. This 
shows that the completely optically thin solution is approached, but 
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Figure 5. Vertical profiles of radiation fields from the disk, for = 3.5. The solid line represents the profiles going through the centre (R = Opc) while the 
other lines (see legend) represent profiles at different radial distance from the plane (R =3, 8, 15 and 24kpc). The different panels show profiles at different 
wavebands, from the B band to the NIR. The profiles in the different wavebands have been scaled as in Fig. [4] The profile overplotted in red in the bottom right 
panel (the K-band case) represents the dustless solution. 



still not reached. The dustless case (the red profile in Fig. [5} indeed 
shows a maximum at z = 0. 

At larger radii the vertical profiles are very flat in all cases, 
because the typical distance to the emitting regions, irrespective of 
whether they are optically thin or optically thick, does not change 
much. 



3.2.1 The Bulge 

Here we discuss the spatial trends of radiation fields for de Vau- 
couleurs bulges (corresponding to a Sersic index n = 4). Devia- 
tions from these trends in bulges with different Sersic indices will 
be disscussed in Sect. [5] Unlike the disk, the RFs of the bulge in 
the plane of the disk (z = 0) and for the optically thick cases (e.g. B 
band) do not show a constant trend with radial distance, even for the 
most optically thick cases (see Fig. [6j. On the contrary, the radial 
profiles show a very steep decrease at small radii. This is because 
of the steep dependence of stellar emissivity on radius, always en- 
suring a monotonic decrease in the ratio of stellar emissivity to dust 
opacity. By the same token, the same trend can be seen in the ver- 
tical profiles at R = (see Fig.|7]l. 

However, going away from the centre, the vertical profiles 
start to qualitatively resemble the behaviour of the disk (a maxi- 
mum in the energy densities at an intermediate distance above the 
plane). This is because the bulge has, like the disk, a layer of stel- 
lar emissivity going above the dust layer, namely an effective ra- 
dius that is larger than the scale-height of the dust disks. Conse- 
quently, there is also a non-monotonic behavior in the radial pro- 
files at larger vertical distances, as in the case of the disk. At larger 
R and z the RF tend to an inverse square law. 



4 VARIATION OF RADIATION FIELDS WITH 
CENTRAL FACE-ON DUST OPACITY t f b 

4.1 The Thin Disk 

The variation of the radiation fields of the thin disk with T f B is ex- 
emplified in Fig. [8] for radial positions in the plane of the disk as 
well as for vertical positions at R = 0, for a fixed wavelength (B- 
band). The first characteristic of the curves is the decrease of the 
energy density of the RF with increasing r^. This is to be expected, 
as photons coming from a wider range of distances can reach any 
given vantage point in the galaxy in the optically thin case than in 
the optically thick case. 

The radial profiles (left panel Fig. [8} also show the expected 
trend with optical thickness, namely the change from a continu- 
ously decreasing energy density of the RF with increasing radius 
for optically thin cases, to, for optically thick cases, a profile which 
is almost flat in the inner regions, only markedly decreasing at high 
increasing radii, (see detailed explanations in Subsect. [3.1.1t . This 
is similar with the trend with wavelength for the same spatial loca- 
tions (see Fig.rrjl, since the increase in wavelength is followed by a 
decrease in the dust opacity at that wavelength. 

The vertical profiles (right panel Fig. [8) show a transition from 
a continuously decreasing energy density of the RF with increasing 
vertical distance from the plane in the optically thin case, to, for 
optically thick cases, profiles which also decrease at smaller dis- 
tances above the plane, but become flat at higher vertical positions. 
This is the expected trend when moving from an optically thin to 
an optically thick case, as explained in Subsect. [3.1.2| (see also the 
similar trend with wavelength in Fig. [3}. 

4.2 The disk 

As in the case of the thin disk, the radial profiles in the plane of 
the disk and the vertical profiles in the centre of the disk (given for 
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Figure 6. Radial profiles of radiation fields from the bulge, for = 3.5. The solid line represents the profiles in the plane (z = Opc) while the other lines 
(see legend) represent profiles at different vertical distances from the plane (z =50, 200, 1000 and 2000 pc). The different panels show profiles at different 
wavebands, from the B-band to the NIR. The profiles have been scaled such that the maximum intensity of the radiation fields of the model galaxy at different 
wavelenths reaches the value of 14.33 X 10 41 erg/pc 3 /A, which is the maximum intensity in the B-band. The calculations are for de Vaucouleurs bulges. 




Figure 7. Vertical profiles of radiation fields from the bulge, for = 3.5. The solid line represents the profiles going through the centre (R = Opc) while the 
other lines (see legend) represent profiles at different radial distance from the plane (R =3, 8, 15 and 24kpc). The different panels show profiles at different 
wavebands, from the B-band to the NIR. The profiles at different wavebands have been scaled as in Fig. [6] The calculations are for de Vaucouleurs bulges. 



the B-band in Fig. [9} show a decrease in the energy density of the 
radiation fields with increasing dust opacity t b . 

The radial profiles for the optically thin cases show a contin- 
uous decrease with increasing radius (left panel Fig. |9), since the 
disk is in these cases optically thin even in the inner regions. In the 
optically thick cases the radial profiles are flat in the inner radii, or 
even slightly increase, followed by a decrease of the energy den- 



sities at larger radii. The disk remain essentially flat for radial dis- 
tances for which the disk is still optically thick. This means that the 
higher the face-on optical depth t b , the larger the radius for which 
the disk is optically thick, and therefore the profiles are flat. 

As discussed in Sec. 13.21 the vertical profiles (right panel 
Fig. [9} show a maximum in the energy density at intermediate ver- 
tical distances, due to the fact that the stellar emissivity in the disk 
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Figure 8. Left: B band radial profiles from the thin disk calculated in the plane (z = Opc). Right: B band vertical profiles from the thin disk calculated through 
the centre (R = Opc). The different lines represent profiles corresponding to different central face-on dust opacities in the B-band (see legend). 
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Figure 9. Left: B 
centre (R = Opc). 



band radial profiles from the disk calculated in the plane (z = Opc). Right: B band vertical profiles from the disk calculated through the 
The different lines represent profiles corresponding to different central face-on dust opacities in the B-band (see legend). 



has a larger scaleheight than the scaleheight of the two dust disks. 
The vertical distance at which this maximum occurs increases with 
increasing dust opacity t£, since the optical thick regime extends 
to higher vertical distances in an optical thick case. This effect is 
particularly strong for the T f B = 8 case, where the energy densi- 
ties increase by a factor of ~ 6, when moving from the centre to a 
position 600 pc above the plane. 

4.3 The Bulge 

The radiation fields of the bulge also decrease with increasing r B , 
albeit at a much smaller rate that for the thin disk and disk (see 
Fig. |T0t . The radial profiles of RFs in the plane of the disk show 
a monotonic decrease with increasing radial distance, irrespective 
of T^j. Similarly, the vertical profiles of RFs at R = Opc show a 
monotonic decrease with increasing vertical distance for all dust 
opacities considered here. As explained in Sect. I3.2.1I even in the 



very optically thick cases the profiles show a steep decrease with R 
and z in the central regions due to the very steep dependence of the 
stellar emissivity with radial and vertical distance. 



5 VARIATION OF RADIATION FIELDS WITH BULGE 
SERSIC INDEX 

In Sects, [3^2Tl and |431 we discussed the spatial variation of the ra- 
diation fields from bulges, as well as the variation with wavelength 
and central face-on dust opacity. This was done for the concrete ex- 
ample of de Vaucouleurs bulges. Here we investigate whether the 
trends seen in the RF for bulges change when different Sersic in- 
dices are considered. In Fig.QT]we show radial and vertical profiles 
of RFs at z = and R = respectively, for four different values 
of Sersic indices, n = 1,2,4, 8, for the case Tg = 3.5, and for the 
B-band. 
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Figure 10. Left: B band radial profiles from the bulge calculated in the plane (z = Opc). Right: B band vertical profiles from the bulge calculated through the 
centre (R = Opc). The different lines represent profiles corresponding to different central face-on dust opacities in the B-band (see legend). The calculations 
are for de Vaucouleurs bulges. 
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Figure 11. Left: B band radial profiles from the bulge calculated in the plane (z = Opc). Right: B band vertical profiles from the bulge calculated through the 
centre (R = Opc). The different lines represent profiles corresponding to different Sersic indices of the stellar emissivity (see legend). The profiles are scaled 
such that the maximum of the energy density of the radiation fields of a model galaxy reaches the value of 39.48 X 10 41 , which is the maximum intensity for 
the model with n = 8. The calculations are for tL = 4.0. 



The radial profiles show a monotonic decrease with increas- 
ing radial distance, irrespective of the Sersic index. However, for 
higher Sersic indices the decrease is steeper, in particular in the 
very central regions, as expected due to the steeper dependence of 
the stellar emissivity. For n = 1 the profiles start to resemble those 
from the disk, due to the exponential form of the stellar emissivity. 

The vertical profiles show more pronounced differences in the 
trends for the different Sersic indices. Thus, for lower Sersic in- 
dex the profiles become less steep. In the extreme case of n = 1 
the profile not only flattens, but even shows a second maximum. 
This exponential form of the stellar emissivity produces a solution 
closer to the disk solution, where a maximum occurs at a certain 
vertical distance from the disk. This is due to the fact that the stel- 
lar emissivity has a larger scale-height than that of the two dust 



disks. However, unlike the disk, there is another maximum in the 
plane of the disk. This is because, although the stellar emissivity is 
described by an exponential distribution, this emissivity is spatially 
distributed in an ellipsoid rather than a disk. 

Overall we can conclude that the radiation fields from bulges 
strongly depend on the assumed Sersic index of the bulge. Because 
of this we include in the library of radiation fields for bulges the 
calculations for all four values of Sersic indices shown above. 



6 THE LIBRARY OF RADIATION FIELDS 

We have created a library of diffuse radiation fields, calculated on 
a cylindrical grid (R, z). The radiation fields are given for different 
values of central face-on dust opacity in the B-band t^, which is 
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the only parameter of the model shaping the spatial variation of the 
RFs. Since the energy densities of the radiation fields are additive 
quantities, they scale with the spatially integrated luminosity den- 
sity at a given wavelength. Therefore we only calculate RFs for a 
fixed reference luminosity density. The radiation fields are sepa- 
rately calculated for the disk, thin disk and bulge. For the bulge we 
considered calculations for different values of the Sersic index of 
the bulge. 

In total we performed calculations for seven values of the cen- 
tral face-on dust opacity, = 0.1,0.3,0.5, 1.0,2.0,4.0,8.0. Solu- 
tions for other values of t^, such as for the case t b = 3.5 used 
to illustrate the solutions in this paper, can be found by interpola- 
tion. We considered four values for the Sersic index of the bulge, 
n = 1, 2, 4, 8. At each wavelength we therefore have 7 (for t b ) x 2 
(for disk and thin disk) + 7 (for t b ) x 4 (for Sersic indices of bulge) 
= 42 combinations at each sampled UV/optical/NIR wavelength. 
We considered calculations at 15 different wavelengths, ranging 
from 912 A to 5 //m, as listed in Table E.2 of Popescu et al. (201 1). 
Thus, the library contains a total of 42 x 15= 630 files with two- 
dimensional spatial grids of energy densities of the RFs. 

As mentioned before, the radiation fields have been calculated 
for a reference spatially integrated luminosity density at each wave- 
length and for each stellar component. The reference luminosity is 
defined according to two main parameters of the model of Popescu 
et al. (2011), the star-formation rate S FR and old. The S FR is a 
parameter defining the luminosity of the young stellar populations 
and old is a parameter defining the luminosity of the old stellar 
population. Because the reference luminosity densities have been 
defined to correspond to the unit values of the two model parame- 
ters, S FR = lM e /yr and old = 1, they are referred to as unit in- 
trinsic spectral luminosity densities of the young and and old stellar 
populations, ^°„™f and L°!f mjl , as defined in Sect. 2.3 and Table E.2 
of Popescu et al. (201 1). Correspondingly, because the energy den- 



Zref 



= Z * 



sities of the radiation fields are calculated for the unit L ' 



and 



, old 



they are referred to as the unit radiation fields, as defined in 



Sect. 2.5.2 of Popescu et al. (2011): U 



disk n ! disk 
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Thus, 
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u ! ., are calculated for a stellar disk having L 

Admit o ! 

calculated for a bulge having L°'f mjr The library of radiation fields 
presented in this paper and made available at CDS is therefore given 
in the form of unit radiation fields. 

The unit radiation fields are calculated for a reference disk 
size, as defined in Popescu et al. (2011). The reference size is 
h.s.ref = 5670 pc, defined as the intrinsic radial scale-length of the 
volume stellar emissivity of the disk in the B-band. For morpholog- 
ical components of galaxies with different luminosities Lf sk ,L t f' sk , 
L b " lge and disk sizes h s , the radiation fields at a position (R,z) can 
be scaled according to the formula 
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When combined with the solutions for the UV/optical/submm 
SEDs found using the library of Popescu et al. (2011), the ampli- 
tude of the radiation fields can be found by scaling the unit radiation 
fields according to the model parameters S FR, old, bulge-to-disk 
ratio B/D and the dumpiness factor F, using the prescription from 
Sect. 2.5.2 in Popescu et al. (201 1). This allows observables calcu- 
lated from the radiation fields, such as gamma-rays or radio emis- 
sion, to be self-consistently combined with the solutions for the 
integrated SEDs, thus expanding the range of applicability of the 
radiation transfer model to a broader range of energies and physi- 
cal quantities. 

As noted in Popescu et al. (201 1), in any dusty galaxy h s will 
not actually be a directly observable quantity, even if the distance 
to the galaxy is known. As described and quantified in Mollenhoff 
et al. (2006) and in Pastrav et al. (2013), the apparent size of a 
galaxian disk obtained from photometric fits to optical images will 
differ from the true size due to the stronger attenuation of light at 
the centre of galaxies compared to the outer regions. Furthermore, 
when bulge-disk decomposition is performed on dusty galaxies, 
there will be additional effects due to dust on the decomposition 
itself (Pastrav et al. in prep), resulting in additional differences be- 
tween the intrinsic and the measured sizes of the disk. Thus, when 
scaling radiation fields, the measured scale-length (or effective ra- 
dius) of the disk should be corrected for the effects of dust. Cor- 
rections for these effects can be found in Pastrav et al. (2013) and 
Pastrav et al. in prep. Since the corrections from Pastrav et al. were 
derived using the same radiative transfer model as used for the cal- 
culation of the radiation fields, incorporation of these corrections 
in the scaling of the RFs allows for a self-consistent analysis of the 
integrated SEDs, structural parameters and radiation fields. 



7 APPLICATION: THE EFFECT OF DIFFERENT BULGE 
SERSIC INDICES ON DUST HEATING 

One important conclusion of this study is that the radiation fields 
from bulges strongly depend on the value of the Sersic index of the 
stellar emissivity (see Sect. [5j. This raises the possibility that the 
heating of dust in the disk may also be strongly influenced by the 
type of bulge a galaxy would host. 

To check this effect we performed calculations for dust heat- 
ing and emission for a model galaxy having bulges described by 
different Sersic indices (n = 1,2,4,8). We chose a typical value 
for the bulge-to-disk ratio, of B/D = 0.33. We also considered a 
galaxy having the unit luminosities (as described in our model - 
see Sect. [6] and Popescu et al. 2011) for the old and young stellar 
populations (S FR = 1 M Q /yr, old = 1) and r B = 4.0. The resulting 
integrated SEDs are shown in the top panel of Fig. [12] One can see 
that with higher the Sersic index, the larger the emission in the 24 
to 60 /im wavelengths range is. This is because of the very strong 
heating of the dust in the very central region of the disk by optical 
photons emitted by the bulge. Although the total dust luminosity 
is almost unaffected by the shift in the Sersic index, the integrated 
mid-infrared/far-infrared colours are somewhat boosted for galax- 
ies containing bulges with high Sersic index. This boost in colour is 
however moderate in comparison with the boost due to an increase 
in the B/D ratio, as shown in Fig. 19 of Popescu et al. (201 1). 

The stronger effect though appears in the surface-brightness 
distribution, in particular in the central regions of galaxies contain- 
ing bulges. In the bottom panel of Fig. [T2] one can see that the 



13 



8 SUMMARY 




10.0 



Figure 12. Spatially integrated diffuse dust emission SEDs (top) and 24 /im 
radial profiles along the major-axis (bottom) of a model galaxy (containing 
a young stellar disk, an old stellar disk - with a B-band exponential scale- 
length of 5.67 kpc - and a bulge with R e ff = 1.3 kpc). The parameters of 
the model were set to SFR = 1 M Q /yr, old = 1, r£ = 4.0, B/D = 0.33. The 
different lines correspond to different Sersic indices of the stellar emissivity 
in the bulge (see legend). The profiles are for a galaxy seen face-on. 



In this paper we provide a library of energy densities of the dif- 
fuse radiation fields arising from stars in the main morphological 
components of spiral galaxies: disks, thin disks and bulges. The 
library contains 630 two-dimensional spatial grids (in cylindrical 
coordinates) of energy densities of the RFs, spanning seven values 
of face-on central optical depth tL, 15 values of wavelength, and, 
for bulge RFs, four values of Sersic index n. Thus, we considered 
T"jg = 0.1,0.3,0.5, 1.0,2.0,4.0,8.0, n = 1,2,4,8 and wavelengths 
from 912 A to 5 fim. The radiation fields are sampled on an irregu- 
lar grid with sampling intervals in the radial direction ranging from 
50 pc in the centre to up to 2 kpc in the outer disk and in the verti- 
cal direction ranging from 50 pc in the plane to up to 500 pc in the 
outer halo. 

The model radiation fields are those used to generate the large 
library of dust- and PAH-emission SEDs presented in Popescu et 
al. (2011). This library is self-consistently calculated with the cor- 
responding library of dust attenuations first presented in Tuffs et 
al. (2004) and then in updated form in Popescu et al. (2011). All 
calculations were made using a modified version of the ray-tracing 
radiative transfer code of Kylafis & Bahcall (1987), and the dust 
model from Weingartner & Draine (2001)and Draine & Li (2007), 
incorporating a mixture of silicates, graphites and PAH molecules. 

We discussed the radial and vertical variation of the RFs, as 
well as their variation with t^, and, for bulges, with the Sersic in- 
dex. We give an analytic formula for the optically thin case of an 
infinitely thin disk, which we use to explain the trends in the radial 
profiles of thin disks. We also show how the library can be self- 
consistently combined with the results from modelling the spatially 
integrated UV/optical/submm SEDs of galaxies from Popescu et al. 
(2011). 

As a practical application we calculated the dust emission 
SEDs for galaxies having different bulge Sersic indices. We find 
strong mid-infrared localised emission in the central regions of 
disks with high Sersic index bulges, which can mimic a central 
star-formation region or a mild AGN activity. Unlike these latter 
two cases, the dust is in this case heated by optical photons from 
the bulge. 



radial profiles of galaxies containing bulges described by different 
Sersic indices can vary by a factor of more than 10 at 24 jim. This is 
remarkable, since dust emission at 24 jim is commonly interpreted 
in terms of heating either by UV photons in star-formation regions 
or from accretion sources (AGN), whereas here the heating is dom- 
inated by old stars in the bulge. Moreover, one can see that this 
enhancement occurs only within a few hundred parsec of the cen- 
tre. At these linear scales, for resolved galaxies beyond 10 Mpc the 
excess emission due to the bulge heating will appear as a sepa- 
rate warm point-like source at the angular resolution of the Spitzer 
Space Telescope at 24 fim. Particular care must therefore be taken 
in interpreting central sources seen in the mid-infrared towards spi- 
ral galaxies with bulges, to distinguish between an obscured low- 
luminosity AGN, star-formation activity and the effect of the bulge 
we have identified here. 
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